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Abstract 

An analytic, parameter-free (up to overall scale factors) solution of the Bohr Hamiltonian 
involving axially symmetric quadrupole and octupole deformations, as well as an inhnite 
well potential, is obtained, after separating variables in a way reminiscent of the Variable 
Moment of Inertia (VMI) concept. Normalized spectra and B{EL) ratios are found to agree 
with experimental data for ^^®Ra and ^^®Th, the nuclei known to lie closest to the border 
between octupole deformation and octupole vibrations in the light actinide region. 


1. Introduction 

Gritical point symmetries Qq are attracting recently considerable interest, since they 
provide parameter-independent (up to overall scale factors) predictions supported by exper¬ 
iment |S1 in El El • The E(5) [T] and X(5) [2] critical point symmetries have been obtained 
from the Bohr Hamiltonian [7] after separating variables in different ways and using an 
inhnite square well potential in the (3 (quadrupole) variable, the latter corresponding to the 
critical point of the transition from quadrupole vibrations [U(5)] to axial quadrupole de¬ 
formation [SU(3)] j2j- A systematic study of phase transitions in nuclear collective models 
has been given in PEiiin]. 

In the present work a solution of the Bohr Hamiltonian aiming at the description of the 
transition from axial octupole deformation to octupole vibrations in the light actinides m 
is worked out. In the spirit of E(5) and X(5) the solution involves an inhnite square well 
potential in the deformation variable and leads to parameter-free (up to overal scale factors) 
predictions for spectra and B{EL) transition rates. Both (axially symmetric) quadrupole 
and octupole deformations are taken into account, in order to describe low-lying negative 
parity states related to octupole deformation, known to occur in the light actinides m 
Separation of variables is achieved in a novel way, reminiscent of the Variable Moment of 
Inertia (VMI) concept ^21- The parameter-free predictions of the model turn out to be 
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in good agreement with experimental data for ^^®Ra ra and jHj, the nuclei known 

jiniiisj to he closest to the transition from octupole deformation to octupole vibrations in 
this region. 

In Section 2 of the present work the solution is worked out. Numerical results are given 
and compared to experimental data in Section 3, while Section 4 contains discussion of the 
present results and plans for further work. 

2. The model 

When only axially symmetric quadrupole {P 2 ) and octupole (/ds) deformations are taken 
into account, the Hamiltonian reads imiiHi 




^ ^^3 2Bx df3x dpx 


"a’ ® 


+ 




6(H2/9| + 253/91) 


+ V{P 2 ,P 3 ) 


( 1 ) 


where B 2 , B^ are the mass parameters. 

One then seeks solutions of the Schrodinger equation of the form im 


<^t{f32,f33,e) = i(52f3s)-^/^^t{f32,f33)\LM0,±), 


( 2 ) 


where M is the angular momentum projection onto the laboratory-hxed 5-axis, = 0 
is the projection onto the body-hxed 5'-axis, and the functions |LM0,-|-) and |LM0, —) 
transform according to the irreducible representations (irreps) A and Bi of the group D 2 
respectively mi CHI, their explicit form being given in [iniiini. 

Introducing miiiH] 
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B 2 + B^ 
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(3) 


reduced energies e = {2B/fi‘^)E and reduced potentials u = {2B/h‘^)V [H Ej, as well as 
polar coordinates (with 0 < /3 < cx) and —tt/2 < 0 < tt/2) [T71ITH] 

^2 =/3 cos 0, p 3 = (3sm(j), p = + Pi, (4) 


the Schrodinger equation takes the form mi 
52 Id L{L + 1) 1 


dp^ pdp 3/?2(l + sin2 0) /?2d02 


+ u{P,p) + 


P"^ sin^ 2p 


— 




(5) 

Separation of variables in Eq. © can be achieved by assuming the potential to be of 
the form [121120] u{P,p) = u{P) +u{p)/P‘^, leading to 

d^ 1 d 


/9^ 


+ I (/^) = 


d /32 p 


( 6 ) 
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m(0) - Ml(0) 1 X^(0) = 


where 


ml(0) = 


+ 


i:(L + i) 


(7) 


( 8 ) 


sin^20 3(l + sin2 0)’ 

with 1/2 being the separation constant and \I/^(0,0) = 0^ (/3)x^(0)) where, however, the ± 
indices have become redundant. 

The potential ul{4>) of Eq. (jH)) is shown in Fig. 1 for several values of L, normalized to 
its minimum value for each L. It is clear that in each case the potential has the form of a 
deep well, possessing an L-dependent minimum, denoted by 0^, and determined from the 
equation 

9 t, -I- 1 1 «in rh 19 pn« Orh 

(9) 


! 2 L(L + 1) sin0cos0 12 cos 20 

'^LWJ ~ yi j • 2 i\2 • 3 r)i ~ 

o (1 + sm cpy sm 20 


Using standard trigonometric identities and dehning x = sin^0 and b = 9/(4L(L + 1)) one 
easily sees that Eq. Q takes the form 


— 2(1 + b)x^ + (1 — 36)a;^ + 6 = 0, 


( 10 ) 


which turns out to have only one real root in the interval 0 < x < 1 (imposed by x = sin^ 0). 

Given its form, the potential ^^(0) can be approximated around the minimum by the 
hrst terms of the Taylor expansion as 


Ml(0) Ml(0l) + ^^[f^^ (0 - 0 l)^. 


( 11 ) 


In Eq. ((7j) one can then omit the potential m( 0), treating ml( 0) as an effective potential 
naturally occuring in the framework of the theory, leading Eq. © into the harmonic 
oscillator form 


92 X 
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+ 


with 


c' = 


/<(0i 


-(0 — 0l)^, El — 


Z /2 - Ml( 0 l) 


u'Ml) 


Since = 2n + 1, where n is the number of oscillator quanta, one obtains 


z/2 = 


lnU<pL 


-{2n + 1 ) + ml(0l). 


( 12 ) 


(13) 


(14) 


In what follows we are going to be limited to the case n = 0. 
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Returning to Eq. using for u{(3) an infinite well potential {u{(3) = 0 if /3 < (3w] 
u{(3) = oo a (3 > Pw), and defining j2] = ^|, ^ = one is led to the Bessel equation 


(Ppu 1 dp; 


+ 


dz‘^ z dz 


+ 


1 - 


u 


pu = 0 . 


(15) 


Then the boundary condition pu{Pw) 


0 determines the spectrum 


^/3,s,u 


= e 


I3,S,L 


= (ks. 
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h — 




Pw 


and the eigenfunctions 

Ps,l'{,P') Ps,l{^P') Cs^i,Jy{ks^i,PP 


( 16 ) 


(17) 


where Xs,y is the sth zero of the Bessel function Jv{z), while Cs^v are normalization constants, 
determined from the condition /(f'^ \ps^y{P)p‘PdP = 1 to be Cg^u = y/2/Jyj^pks^v)- The 
notation has been kept similar to that of Ref. |2] . 

A few comments are now in place: 

a) L-dependent potentials, as the one of Eq. (jSl), are known to occur in nuclear physics in 

the framework of the optical model potential EH. as well as in the study of quasimolecular 
resonances, such as j 22 ]. 

b) The procedure followed for the determination of pi is reminiscent of the Variable 
Moment of Inertia (VMI) model ^2] • In the VMI case the energy is minimized with respect 
to the moment of inertia for each L separately, resulting in a moment of inertia increasing 
with L. In the present case the effective potential energy ul is minimized with respect to 
0 for each L separately, resulting in pi values increasing with L. As a consequence, in Eq. 
(jH|l the denominator of L{L + 1), which can be considered roughly as playing the role of a 
moment of inertia, is also increasing with L. 

c) In order to separate variables, one can in general assume u^P,^) = u{P) +u{'y), as 
in the X(5) model j2], or u{P,'y) = u{P) + u{'y)/P'^, as in Refs. [T71 OH] . In the former case 
the separation is approximate, since a term is involved in the 7 equation, replaced by its 
average value, but no extra parameter is introduced in the /3-equation by the separation. 
In the latter case, the separation is exact, but (at least) one extra parameter appears 
in the /3-equation, coming in from the 7 -equation through the separation constant. In the 
present model this disadvantage of the latter case is avoided through the VMI-like procedure 
adapted. 

d) For each L the specific value of the variable p, which decides the relative presence 
of the quadrupole and octupole deformations, is determined in Eq. m by the effective 
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potential ul{4)), which has a rigid shape, as seen in Fig. 1, while the potential m( 0) plays 
no role. 

The calculation of B{EL) transition rates proceeds as in Ref. and need not be 
repeated here. 

3. Numerical results and comparison to experiment 

The parameter-free predictions of the model for the lowest bands are given in Table 1, 
together with the experimental spectra of ^^®Ra ^ 2 ] and ^^®Th ^ 3 ], which are known to 
lie near the border between the regions of octupole deformation and octupole vibrations 
uniiiE], as also seen in Figs. 2(a) and 2(b), where the experimental spectra of ^^®Ra 
and 220-234rpj^ included. In both hgures the region below the theoretical predictions cor¬ 
responds to octupole deformation, characterized by minimal odd-even staggering rapidly 
decreasing and disappearing with increasing L, while the region above the theoretical pre¬ 
dictions corresponds to octupole vibrations, characterized by large odd-even staggering 
decreasing very slowly with increasing L. As seen in Table 1, in the case of ^^®Ra and ^^®Th 
the odd-even staggering is non-negligible only for the lowest four odd levels, the agreement 
between theory and experiment being very good for the even levels, as well as with the rest 
of the odd ones. The absence of staggering in the present model is due to the fact that the 
(inhnite) potential wells for /ds < 0 and /ds > 0 [see Eq. (jH)] are separated by an inhnite 
barrier, and not by a hnite one, as needed for odd-even staggering to be present Pi. 

Several parameter-free predictions for B{E1), B{E2), and B{E3) transition rates, ap¬ 
propriately normalized, are reported in Table 2 . Since the lack of experimental data does 
not allow for direct comparison to experiment, comparisons to B{E1) branching ratios of 
^^®Ra |22] and B{E1)/B{E2) ratios for ^^®Th jT^l are shown in Figs. 3(a) and 3(b) respec¬ 
tively. The theoretical predictions in Fig. 3(a) are compatible with the data within the 
experimental errors, while in addition they are quite similar to predictions of the Extended 
Coherent States Model (ECSM) [221 obtained with the lowest order El operator (R-h), as 
well as with two different choices of the El operator including anharmonicities (R-I, R-II) 
m- The theoretical predictions in Fig. 3(b) are also compatible with the data within the 
experimental errors, lying considerably lower than the predictions of Ref. Hg (BBS). 

It is worth remarking that the ground state band spectrum and intraband B{E2)s of 
the present model are quite similar to these of the X(5) model [2j (extensively tabulated in 
Ref. |2Hj)- Indeed the present model can be considered as an extension of X(5), in which 
the octupole degree of freedom is taken into account in order to account for the low-lying 
negative parity bands, while in parallel the 7 degree of freedom is left out in order to keep 
the problem tractable. One important difference between the two models is related to the 
(normalized) position of the 0^ state, which is predicted at 5.649 by the X(5) model j2]. 
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but at 12.569 by the present model. This implies that while searching for X(5)-like nuclei 
in the light actinide region, one should expect the 0^ state to occur higher by a factor of 
two. Indeed the Ra and Th isotopes near A = 226 exhibit high-lying 0^ states nsuni. 

It should also be noticed that the predictions of the present parameter-independent 
model are very similar to those of the one-parameter ((/)o) Analytic Quadrupole Octupole 
Axially symmetric (AQOA) model in which best agreement to experiment is obtained 
for 00 = 56° in the case of ^^®Ra, and for 0o = 60° in the case of ^^®Th. These values of 0o 
are understandable, when compared to the 01, values shown in Table 1. 

4. Discussion 

A parameter-free (up to overall scale factors) version of an Analytic Quadrupole Oc¬ 
tupole Axially symmetric (AQOA) model involving an inhnite well potential, suitable for 
describing the transition from octupole deformation to octupole vibrations in the light ac¬ 
tinides, has been constructed, after separating variables in the Bohr Hamiltonian in a way 
reminiscent of the Variable Moment of Inertia concept. Spectra and B{EL) ratios are shown 
to be in good agreement with experimental data for ^^®Ra and ^^®Th, the nuclei supposed 
to lie closest to the above mentioned border. Application of the model in the A k, 150 
region, where octupole deformation is also well established nn, is receiving attention. 
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Table 1: Spectra of the present model for the ground state band and the associated negative 
parity band (s = 1), as well as for the hrst excited band (s = 2), together with relevant val¬ 
ues of (pL and experimental data for ^^®Ra and ^^®Th m. Each spectrum is normalized 
to the energy of its own 2^ state. See section 3 for further discussion. 



0 L 

th 

^^«Ra 

226Th 


0 L 

th 

^"«Ra 

226Th 

s = 1 
0 + 

45.00 

0.000 

0.000 

0.000 

1 - 

45.70 

0.337 

3.747 

3.191 

2 + 

47.01 

1.000 

1.000 

1.000 

3“ 

48.77 

1.967 

4.749 

4.259 

4+ 

50.77 

3.200 

3.127 

3.136 

5- 

52.79 

4.657 

6.603 

6.240 

6 + 

54.71 

6.297 

6.155 

6.195 

7- 

56.47 

8.093 

9.264 

9.112 

8 + 

58.05 

10.025 

9.891 

9.999 

9- 

59.46 

12.081 

12.677 

12.785 

10 + 

60.73 

14.254 

14.185 

14.409 

11 - 

61.86 

16.537 

16.743 

17.152 

12 + 

62.89 

18.928 

18.931 

19.324 

13- 

63.81 

21.424 

21.388 

22.105 

14+ 

64.65 

24.023 

24.061 

24.675 

15- 

65.42 

26.724 

26.536 

27.554 

16+ 

66.13 

29.526 

29.523 

30.413 

17- 

66.78 

32.427 

32.126 

33.418 

18+ 

67.38 

35.428 

35.300 

36.497 

19- 

67.94 

38.527 

38.099 

39.627 

20 + 

68.46 

41.724 

41.375 

42.896 






s = 2 
0 + 

45.00 

12.569 

12.186 

11.152 






2 + 

47.01 

14.253 








4+ 

50.77 

17.871 
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Table 2: B{EL\ Li Lf) values between states of the present model with s = 1. B{E2)s 
with Li and Lf even are normalized to the 2f 0^^ transition, while B{E2)s with Li and 
Lf odd are normalized to the Sf —If transition. B{El)s are normalized to the Ij" —0^" 
transition, while B{E3)s are normalized to the 3^ —Oi" transition. See section 3 for further 
discussion. 


LI 

rn 

B{E2) 

LI 

r TT 

B{E1) 

LI 

r TT 

B(E3) L- 

r TT 

B{E3) 


2 + 

0 + 

1.000 

1 - 

0 + 

1.000 

3- 

0 + 

1.000 

2 + 

1 - 

1.793 

4+ 

2 + 

1.494 

2 + 

1 - 

1.238 

4+ 

1 - 

1.358 

3- 

2 + 

1.362 

6 + 

4+ 

1.749 

3- 

2 + 

1.389 

5- 

2 + 

1.581 

4+ 

3- 

1.341 

8 + 

6 + 

1.944 

4+ 

3- 

1.522 

6 + 

3- 

1.752 

5- 

4+ 

1.369 

10 + 

8 + 

2.103 

5- 

4+ 

1.650 

7- 

4+ 

1.896 

6 + 

5- 

1.411 

12 + 

10 + 

2.237 

6 + 

5- 

1.775 

8 + 

5- 

2.022 

7- 

6 + 

1.456 

14+ 

12 + 

2.350 

7- 

6 + 

1.896 

9- 

6 + 

2.136 

8 + 

7- 

1.500 

16+ 

14+ 

2.448 

8 + 

7- 

2.012 

10 + 

7- 

2.238 

9- 

8 + 

1.542 

18+ 

16+ 

2.534 

9- 

8 + 

2.123 

11 - 

8 + 

2.331 

10 + 

9- 

1.582 

20 + 

18+ 

2.611 

10 + 

9- 

2.228 

12 + 

9- 

2.416 

11 - 

10 + 

1.620 




11 - 

10 + 

2.328 

13- 

10 + 

2.494 

12 + 

11 - 

1.655 

3- 

1 - 

1.000 

12 + 

11 - 

2.424 

14+ 

11 - 

2.566 

13- 

12 + 

1.689 

5- 

3- 

1.246 

13- 

12 + 

2.515 

15- 

12 + 

2.632 

14+ 

13- 

1.720 

7- 

5- 

1.413 

14+ 

13- 

2.602 

16+ 

13- 

2.694 

15- 

14+ 

1.750 

9- 

7- 

1.547 

15- 

14+ 

2.685 

17- 

14+ 

2.752 

16+ 

15- 

1.777 

11 - 

9- 

1.658 

16+ 

15- 

2.764 

18+ 

15- 

2.806 

17- 

16+ 

1.804 

13- 

11 - 

1.752 

17- 

16+ 

2.840 

19- 

16+ 

2.856 

18+ 

17- 

1.829 

15- 

13- 

1.832 

18+ 

17- 

2.913 

20 + 

17- 

2.904 

19- 

18+ 

1.852 

17- 

15- 

1.902 

19- 

18+ 

2.983 




20 + 

19- 

1.875 

19- 

17- 

1.964 

20 + 

19- 

3.050 
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Figure 1: Potential uz, [Eq. (jH))] for different values of the angular momentum L, normalized 
for each L to its minimum value. See section 2 for further discussion. 




angular momentum L angular momentum L 


Figure 2: (a) Experimental energy ratios R{L) = E{L)/E{2i) for ^^^Ra [3122], 

[3, and |T2|, compared to theoretical predictions, (b) Same for 220-228rpp 

230Th [3, 232Th [313, and 234Th 3. 
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Figure 3: (a) Experimental B{E1]L L + 1) / B{E1]L —> L — 1) ratios [23] of B{E1) 
values originating from the same level of ^^®Ra, compared to three different theoretical 
predictions from Ref. |2Zj (labeled as R-h, R-I, R-II), as well as to theoretical predictions of 
the present model. See section 3 for further discussion, (b) Experimental B{E1] L ^ L — 1) 
/ B{E2; L ^ L — 2) ratios (multiplied by 10^) |T3| of B[E1) and B(E2) values originating 
from the same level of ^^®Th, compared to theoretical predictions by Bizzeti and Bizzeti- 
Sona [T3] (labeled as BBS), as well as to predictions of the present model. The ratios 
corresponding to L = 10 and 11 have been used for normalization, as in Ref. ra¬ 
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